The entropy puzzle and the quark combination model 
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We use two available methods, the Duhem-Gibbs relation and the entropy formula in terms 
of particle phase space distributions, to calculate the entropy in a quark combination model. The 
entropy of the system extracted from the Duhem-Gibbs relation is found to increase in hadronization 
if the average temperature of the hadronic phase is lower than that of the quark phase. The increase 
of the entropy can also be confirmed from the entropy formula if the volume of the hadronic phase 
is larger than 2.5-3.0 times that of the quark phase. So whether the entropy increases or decreases 
during combination depends on the temperature before and after combination and on how much 
expansion the system undergoes during combination. The current study provides an example to 
shed light on the entropy issue in the quark combination model. 



The quark combination model [U, I2I] (QCM) or its 
variants such as the quark recombination or coalescence 
model are used to describe multi-hadron processes in high 
energy electron-positron, hadron-hadron and nucleus- 
nucleus collisions. Recently they are successful in ex- 
plaining the data at Relativistic Heavy Ion Collider 
(RHIC), which showed that the combination is more 
preferable than the fragmentation for the hadronization 
mechanism in intermediate transverse momentum region 
in heavy ion collisions at high energies [3-5]. The ba- 
sic idea of the QCM is that hadrons are combined from 
quarks and anti-quarks, i.e. two quarks merge into a me- 
son (a 2-to-l process) and three quarks into a baryon (a 
3-to-l process). One of the most frequently asked ques- 
tions about the validity of the QCM is: does the entropy 
decrease in such 2-to-l and 3-to-l processes? The ques- 
tion arises naturally from the common sense that the en- 
tropy is normally reduced in a particle number decreas- 
ing process. The entropy problem is essential to nucleus- 
nucleus collisions since the issue of the entropy produc- 
tion is closely related to thermalization and has been the 
concern of the researchers for many years. There have 
been previous attempts or arguments about the entropy 
in the QCM- like models in literature but no sys- 

tematic investigation has yet been made. In this paper 
we will address the entropy problem in the course of the 
quark combination in heavy ion collisions. As an exam- 
ple we will use a specific version of the QCM which is 
developed by the Shandong group [l^-O- The advan- 
tage of the Shandong QCM is its exclusive nature and 
simplicity, i.e. it can give all hadrons in an event. Con- 
sidering that the entropy can be regarded as a measure 
of degrees of freedom, the exclusive nature of Shandong 
QCM is suitable and applicable to the entropy problem. 

The building blocks for the QCM are constituent 
quarks and anti-quarks of the light flavors u,d and s with 
the momentum distributions /i(p) = fi{y,PT) jT^ - flGj . 
All hadrons come out with fixed momenta from combi- 
nation of their constituent quarks. Note that our QCM 
(or any of QCM-like models) does not deal with how the 



constituent quarks acquire masses which is closely related 
to the physics of chiral phase transition, although the 
QCM-like models can be regarded as a phenomenologi- 
cal and microscopic description of the phase transition. 
There are no gluons in the model because the role of 
a gluon can be replaced by a pair of quarks. Actually 
hadronization is a dynamic process which takes place lo- 
cally at a freeze-out temperature, while our QCM (or 
any of QCM-like models) lacks the real time feature but 
treats hadronization in a global or an effective way. In 
this sense the QCM is similar to the statistical models. 

Suppose we have Ni quarks or anti-quarks with the fla- 
vor i — u,d,s,u,d, s, the total number of quarks and anti- 
quarks are then Ngg = Nq + Nq with Ng = Nu + Nd + Ns 
and Nq = + N-^ + N-g. In mid-rapidity and high 
energies almost all quarks and anti-quarks are excited 
from vacuum in pairs, we have A^^ w Nq. In lower en- 
ergies there are remenants of nucleons participating in 
the collisions which are left in the mid-rapidity region 
and leads to Nq > Nq. Since the strange quarks are 
heavier they are less produced and we can use a fac- 
tor As — Ng/Nu ~ Ng/Nd < 1 to characterize the 
suppression of strange quarks relative to up and down 
quarks. The procedure to randomly combine all quarks 
and anti-quarks into hadrons is quite straightforward: we 
line up all these quarks and anti-quarks randomly into 
a 1-dimensional queue and let them combine into the 
groud state hadrons following the combination rule [13- 
[l3|. The momentum distribution of a specific hadron is 
the convolution of those of its constituent quarks. Here 
the groud state hadrons include (3/2)+ decuplet and 
(1/2)+ octet baryons and 0~ pseudoscalar and 1~ vector 
mesons. The short-life resonances are allowed to decay 
into the long-life or final state hadrons which are recorded 
in the detectors. The experimental data are used to fix 
the inputs of the model. The charged multiplicity can fix 
Nq and Nq. The ratio of multiplicity of kaons to pions 
can fix Ag. Note that we focus on the unit rapidity re- 
gion at the central rapidity \y\ < 0.5. One of the most 
important element of the model is the transverse momen- 
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turn spectra of thermal quarks which are parameterized 
by the blast wave model ^17i] , 
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where ^ — r/Rmax with Rmax the maximum radius of 
the thermal source, Tf = 165 MeV is the local tempera- 
ture, p = tanh~"'^(/3s^") with jig the radial velocity on the 
surface of the firsball and n = 0.3 a parameter, Jq and 
Ki are the modified Bessel functions. The average radial 
velocity is given by {Pr) = "^^Ps/in + 2). By fitting the 
data for transverse momentum spectra of hadrons we can 
fix the parameters /3s or equivalently (/?,.), see Tab. [Ill 
With these inputs the available data for the transverse 
momentum spectra and rapidity density of hadrons have 
been well described [l3 - [l6j . The reason we only consider 
thermal quarks is that the entropy is a bulk property 
which is dominated by the small px region. 

Equipped by the above settings, we are able to calcu- 
late the entropy of the system before and after hadroniza- 
tion. Given all particle momentum distributions, one way 
to obtain the entropy of a system is through the Duhem- 
Gibbs relation \\m. 
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where S, E, P, V and T are the entropy, energy, pressure, 
volume, and temperature of the system respectively, and 
Pi and Ni are the chemical potential and number for the 
quark/hadron species i respectively. Note that the sys- 
tem will undergo expansion during hadronization which 
is indeed a dynamic and real time process and takes place 
in local volume. The global volume of an expanding sys- 
tem is not well defined in a Lorentz invariant way. Here 
the volume V is only an effective parameter to character- 
ize the average particle occupancy in spatial dimension. 
The chemical potential pi can be determined from the 
multiplicity ratio of the particle i and to its anti-particle 
I via N,/Nj = e^^'/^ where T is set to 170 MeV for aU 
particles. The quantities on the right hand side of Eq. 
(12) are all known with. 
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where p = |p| and — \/ p^ + mf denotes the scalar mo- 
mentum and energy respectively. For constituent quarks 
we choose rriu = rud = 330 MeV and = 550 MeV. 

In Tab. |T] we hst the total E, PV, pN = J2ilJ'i^i 
and TS per unit rapidity in the central rapidity region 



\y\ < 1 for the quark and hadronic phases. We can see 
that {PV)h (for hadrons) is slightly smaller than {PV)q 
(for quarks) because 
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where in the last line is due to that quarks with the same 
momentum direction combine. We also see that pN is 
very small compared to PV and E since in high energy 
collisions almost all quarks in central rapidity are ex- 
cited in vacuum. Note that the energy per unit rapidity 
are not exactly equal for the two phases and the differ- 
ence is within a few percent. The effect of the energy 
non-conservation is the result of the momentum conser- 
vation imposed during the course of a specific combina- 
tion, because energy and mometum conservation cannot 
be simultaneously fuUfilled with fixed quark masses. For 
example two light constituent quarks with the total mass 
660 MeV combine into a pion of mass 140 MeV, the en- 
ergy of pion is less than that of two quarks for small 
momenta. But the non-conservation effect is small on 
average. With {PV)q w {PV)h and Vq < Vh we obtain 
Pq > Pfi and then Pq — B = Ph {B is the bag constant), 
which is consistent to the condition for phase transition. 
Finally we see that TS per unit rapidity is almost con- 
stant for two phases at each collisional energy upto a few 
percent. If we assume a sudden hadronization the tem- 
perature should be the same for the two phases, then we 
have the approimate conservation of the total entropy in 
hadronization. Actually the hadronization time is finite 
during which the system expands and cools down. There- 
fore the total entropy S of the hadronic phase (with lower 
temperature) is larger than that of the quark phase (with 
higher temperature). All above observations are about 
the directly produced hadrons right after the combina- 
tion. For the final state hadrons after resonance decays 
we also see that E, PV, pN and TS do not change much 
as compared to the directly produced hadrons. This ob- 
servation that the resonance decays cannot compensate 
much to the total entropy is different from previous ar- 
guments that resonant decays might be important. 

Another way of obtaining the entropy is directly 
through the formula. 
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where di are degeneracy factors for quark or hadron 
species i, the signs +/— are for bosons/ fermions, 
and gi(r, p) are phase-space distributions satisfying 



3 



Table I; Total E, PV , fiN = J2i and TS per unit ra- 

pidity in the central rapidity region |y| < 1 for the quark and 
hadronic phases. The unit is GeV. The direct hadrons are 
those directly from combination without any resonance de- 
cays. The final hadrons are long-life hadrons including the 
contributions from resonance decays, namely tt*'", , , 
"K**, p, p, n and n. 
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Figure 1: (Color online) The entropy ratio R — SyJSq as 
functions of the volume ratio x — Vh/Vq at four collisional 
energies. 
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Table II: Radial flow parameters in the transverse momentum 
spectra for quarks at four collisional energies. The unit for 
energies is GeV. 



y/SNN 


(Pr) for u,d 


Wr) for s 


17.3 


0.38 


0.416 
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0.483 


0.5 
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0.55 


200 


0.562 
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/ (27r)i' '^idii^i P) ~ fiip)- the classlcal limit where 
(7i(r, p) ^ 1, one reaches, 

/ -7^:33- ["5»(r,P) In ff,:(r,p) +5»(i-,P)] ■ 
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We can take pions which are most populated particles 
as an example to estimate the order of magnitude of 
the phase space distribution. The multiplicity rapid- 
ity density of 7r+ in most central collisions (centrality 
0-5%) at 200 GeV A is about 320. Suppose most pi- 
ons are limited inside a momentum volume of the size 
(Ap)^ 27rp|cosh2/(ApTA2/) - 10 GeV^ and inside 
the spatial volume of about Vu 2200 fm'^ , then we get 
(77r(r, p) ^ 0.027. This provides an upper limit, for other 
particles gi{r,p) <C 1, therefore Eq. ([H]) is a quite good 
approximation of Eq. ([5]) . 

In order to calculate the entropy via Eq. (tSHH]) from the 
known momentum distributions /i(p) in the QCM, we 
make an approximation or give an estimate for (?i{r, p) by 
g,(r,p) « [i2TT)y{d,V)]Mp), where V = Vg,Vh are the 
effective volumes of quark and hadronic phases respec- 
tively. We can then obtain the entropy ratio R = Sh/Sq 
as functions of volume ratio x = Vh/Vq- The volumes 



Vh of the hadronic phase is taken to be the chemical 
freezeout values of the fireball which are extracted from 
the thermal model (l9l . [20| at four collisional energies, 
T4(fm3) =1125 (17.3 GeV), 1620 (62.4 GeV), 1900 (130 
GeV), 2200 (200 GeV). The results are shown in Fig. [H 
For the above set of values of Vh, only if x > 2.5 — 3.0, 
which corresponds to about 1.35-1.44 for the ratio of the 
fireball radii in the hadronic to quark phase, should the 
entropy of the hadronic phase be larger than that of the 
quark one. 

As a contrast we also consider an ideal case for con- 
stituent quarks whose phase space distribution follows 
Fermi-Dirac distribution, gi(r,p) = l/[cxp(ei/T) + 1]. 
The hadron spectra can be determined by combination 
of quarks. Then we use the Duhem-Gibbs relation ([2]) 
and entropy formula ([5)) to obtain the entropies for quark 
and hadronic phases. The result from the Duhem-Gibbs 
relation is {E,PV,TS) =(0.61,0.145,0.755) (quark), 
(0.59,0.136,0.726) (hadron) for unit volume of quark 
phase. We can also see that TS is almost a constant, 
same as in the real case. The result from the entropy for- 
mula shows the same behavior as in Fig. [T] except that 
the entropy of the hadronic phase is larger than that of 
the quark one happens when x > 16. 

A few comments about our approach and results are 
needed. Our results show that the increase of the en- 
tropy requires an effective and adequate expansion of 
the fireball during hadronization, which also implies that 
hadronization takes finite time to complete. Actually 
in hydrodynamic simulation of the evolution of the fire- 
ball, hadronization is indeed a dynamic process with non- 
trivial space-time profile, i.e. it does not take place uni- 
formly in the same space-time but on the freezeout hy- 
persurface, so different part of the fireball hadronizes in 
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different time. Our current approacfi does not take sucli 
a space-time picture but just provide an averaged ef- 
fect which is more simple and transparent than a real 
hydrodynamic simulation. If combination is treated lo- 
cally one can still compute the entropies before and after 
hadronization for each space-time cell using Eq. ([2]) or 
Eqs. (OS]). The total entropy is a sum over all cells and 
the result is similar to the ideal case in our paper. This 
implies that the entropy can be described in a global and 
effective way as in our current model. In other words 
the entropy is a global quantity which should be insensi- 
tive to the local fine structure. As we have emphasized 
in the beginning that we do not address in our QCM 
the entropy issue in the context of phase transition. We 
just made a few comments about it. If local equilirium is 
reached for a closed system the entropy would not change 
during the transition from the quark to hadronic phase, 
it is the entropy density that changes (decreases) during 
the transition acompanied by the volume expansion. If 
the system is not in local equilibrium the phase transi- 
tion is not well defined (it is indeed a crossover) but still 
one can obtain the total entropy which should increase. 
Such a study in the QCM is independent of whether the 
entropy increases or decreases beyond the combination 
process. 

In summary, we have investigated the issue of the en- 
tropy in the framework of the quark combination-like 
model. As an example for such types of models we used 
the one developed by Shandong group whose exclusive 
nature makes a transparent calculation feasible. We used 



two available methods to calculate the entropies for the 
quark and hadronic phases, one from the Duhem-Gibbs 
relation, another from the entropy formula in terms of 
particle phase space distributions. We found that the 
total entropy from the Duhem-Gibbs relation always in- 
creases in hadronization if the average temperature of 
the hadronic phase is lower than that of the quark phase. 
The increase of the entropy during hadronization can also 
be confirmed from the entropy formula if the volume of 
the hadronic phase is larger than 2.5-3.0 times that of 
the quark phase. This implies that the expansion of the 
fireball takes place during hadronization and it takes fi- 
nite time for the quark phase to hadronize. So whether 
the entropy increases or decreases during combination de- 
pends on the temperature before and after combination 
and on how much expansion the system undergoes during 
combination. The current study provides an example to 
shed light on the entropy issue in the quark combination 
model. 
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